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\Q (54) Title: HIGHLY TIME RESOLVED IMPEDANCE SPECTROSCOPY 

(57) Abstract: A highly time resolved impedance spectroscopy that enhances the measurement of the dynamics of non-stationary 
systems with enhanced time resolution. The highly time resolved impedance spectroscopy includes an optimized, frequency rich 
a.c, or transient, voltage signal is used as the perturbation signal, non-stationary time to frequency transformation algorithms are 
used when processing the measured time signals of the voltage and current to determine impedance spectra which are localized 

^ in time; and the system-charachterizing quantities are determined from the impedance spectra using equivalent circuit fitting in a 

^* time-resolution-optimized form. 
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HIGHLY TIME RESOLVED IMPEDANCE SPECTROSCOPY 

RELATED APPLICATIONS 

5 This Application claims priority to pending German patent application 

19949107.0 entitled HIGHLY TIME RESOLVED IMPEDANCE SPECTROSCOPY 
filed October 12, 1999, and pending German patent Application 10018745.5 entitled 
HIGHLY TIME RESOLVED IMPEDANCE SPECTROSCOPY filed April 15, 2000 
both of which are incorporated by reference in their entirety herein. 

10 

BACKGROUND 

1. Field of the Invention 

15 This invention relates generally to the measurement of electrical impedance 

and, more particularly, to electrochemical impedance spectroscopy. 



2. Description of Related Art 



20 Impedance spectroscopy is a procedure used to characterize the electrical and 

electrochemical properties of investigated systems, and their changes over time. 
Typically, an a.c. voltage signal is applied between a working electrode and a counter 
electrode. If applicable, a simultaneously applied d.c. bias voltage is monitored with 
a reference electrode. Both the applied a.c. voltage signal, as well as the current 

25 response of the system, are measured. The complex electrical resistance (termed the 
impedance Z(co)) of a system can be calculated as a function of the frequency from the 
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quotients of the voltage and current signals in the frequency domain according to 
equation (1). The impedance values for a number of frequencies define the 
impedance spectrum. 



/ (to) 



Various electrical properties of the system or electrochemical processes can be 
derived from the characteristics of the impedance spectra. Particularly for systems in 
which direct current cannot flow, a.c. or transient voltage signals must be used for 
investigations. Due to the high information content of impedance spectroscopy, it is 

10 frequently the preferred technique for measurement of impedance spectra. For 
example, in electrochemistry, impedance spectroscopy is a standard analysis 
technique for investigating e.g. corrosion processes, redox reactions, liquid and solid 
electrolytes, thin polymer films, membranes and batteries. Several papers have 
provided an introduction and overview of the technique and application of 

15 electrochemical impedance spectroscopy. See J.R. MacDonald: "Impedance 

Spectroscopy." (John Wiley & Sons, New York: 1987) and C. Gabrielli: Technical 
Report No. No. 004/83. 1983; C. Gabrielli: Technical Report No. part. No. 12860013. 
1990.), both of which are incorporated by reference in their entirety herein. 

Impedance spectroscopy is also used to characterize semiconductor materials 

20 See A. Bard: Electrochemical Methods . (Wiley & Sons, New York: 1980); and in 

biotechnology (See B.A. Cornell, Braach-Maksvytis, L.G. King et al.: "A Biosensor 
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that Uses Ion-Channel Switches." Nature . 387, p.580-583 (1997). S. Gritsch, P. 
Nollert, F. Jahnig et aL: "Impedance Spectroscopy of Porin and Gramicidin Pores 
Reconstituted into Supported Lipid Bilayers on Indium-Tin-Oxide Electrodes." 
Langmuir . 14 (11), 3118-3125 (1998). C. Steinem, A. Janshoff and M. Siber: 
5 "Impedance Analysis of Ion Transport Through Gramicidin Channels Incorporated by 
reference in Solid Supported Lipid Bilayers." Bioelectrochemistry and Bioenergetics . 
42 (2), 213 (1997). All of the above referenced publications are incorporated by 
reference in their entirety herein. 

The use of impedance spectroscopy has increased greatly, particularly in the 

10 field of biotechnology. In most cases, the electrodes are modified by chemical or 
physical coupling of biofunctional molecules and aggregates (e.g. lipid/protein 
membranes). Impedance spectroscopy is also used to detect adsorption processes. 

There are two forms of impedance spectroscopy: Measuring impedance 
spectra in the frequency domain, Method I; and in the time domain, Method II. 

15 Method I (frequency domain procedure): In the first form, a sinusoidal signal 

at a constant frequency and amplitude is applied within a discrete period, and the 
complex impedance of this discrete frequency is determined. To obtain a spectrum, 
sequential signals at different frequencies are applied. The time resolution, defined as 
the length of time in which the determined spectra follow each other, is low in this 

20 form of impedance spectroscopy. The time for acquiring the data records that 
compose the spectrum is a multiple of the period of the lowest frequency contained in 
the spectrum. The precise duration also depends on the number of the frequencies in 
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the spectrum. Following a frequency change, a transition period is allowed for the 
system to attain an equilibrium. The time resolution of a typical sequence of spectra 
is a few seconds to minutes depending on the observed frequency band. 

Method II (time domain procedure): In the second form, a frequency rich a.c. 
5 voltage signal is applied such as square wave pulses, structured or white noise. By 
using Fourier transformation, the impedance spectrum can be determined from a 
single data record of the time course of the voltage and current signal. Therefore, the 
impedance spectrum is limited regarding the bandwidth and frequency resolution by 
the known limitations of Fourier transformation. The measurement time is normally 

10 at least as long as one period of the lowest frequency in the spectrum of interest. 
Usually a measuring period of several periods of the lowest frequency in the spectrum 
is required to sufficiently improve the signal to noise ratio. The maximum time 
resolution depends on the repetition rate at which the data records, or sets, for Fourier 
transformation are acquired. Because the impedance of all frequencies of the 

1 5 spectrum are measured simultaneously in this method, the time resolution is usually 
much better than that of the first method. 

Method I is normally used to characterize stationary systems or systems 
exhibiting slow dynamics. Commercial devices (frequency response analyzers 
(FRA)), are available for these measurements. At present, method II is primarily used 

20 for measurements where the impedance spectra contain very low frequencies, for 
example, down to about 10" 4 Hz, as required in corrosion studies. 
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The electrical properties of non-stationary systems, which means systems 
whose properties are not constant over time, cannot be measured in many cases with a 
sufficient time resolution by either Method I or Method II procedures of impedance 
spectroscopy. The time averaging effect of method I (summing for several periods of 
5 all the frequencies in the spectrum) and method II (over several periods of the lowest 
frequency contained in the spectrum) does not allow changes in the system over time, 
which are faster than the averaging time, to be resolved by a sequence of impedance 
spectra. The averaging time must be greatly reduced for impedance spectrometers to 
measure non-stationary systems with sufficient time resolution. A single impedance 
1 0 spectrum would then indicate the electrical states of the system localized in time. In 
addition, the individual spectra must be determined with high repetition rate to 
determine the time course of the system-characterizing quantities with a maximum 
time resolution. 

An example of a non-stationary system that have not been able to be measured 
15 with conventional impedance spectroscopy includes lipid bilayer membranes with 
integrated, switching ion channels. The kinetics of many biological processes such as 
opening and closing ion channels in lipid bilayer membranes occurs on a time scale of 
a few milliseconds. These systems are highly relevant in the fields of biotechnology 
and human physiology. 

20 Another example of a non-stationary system that has not been able to be 

measured with conventional impedance spectroscopy is metal and semiconductor 
interfaces with liquid and solid electrolytes with highly dynamic interface processes. 
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In characterizing semiconductors and in the field of electrochemistry, conventional 
impedance spectroscopy cannot be used for many dynamic processes such as the in 
situ observation of rapid etching processes or the relaxation of electrochemical 
systems after voltage jumps since the necessary time resolution is impossible in the 
5 required bandwidths. 

From the discussion above, it should be apparent that there is a need for a 
impedance spectroscopy method and apparatus that can measure non- stationary 
systems with high dynamics. The present invention fulfills this need. 



1 0 SUMMARY OF THE INVENTION 

A method and apparatus to measure impedance and impedance spectra in 
rapid sequence. Measurement of impedance, and impedance spectra, in rapid 
sequence provides repeated or continuous characterization of the electrical properties 
of the system under investigation. Highly time resolved impedance spectroscopy 

15 enhances the measurement of the dynamics of non-stationary systems due to its 
enhanced time resolution. 

Three aspects of highly time resolved impedance spectroscopy are: (1) an 
optimized, frequency rich a.c, or transient, voltage signal is used as the perturbation 
signal; (2) non-stationary time to frequency transformation algorithms are used when 

20 processing the measured time signals of the voltage U(t) and current I(t) to determine 
a sequence of impedance spectra, where each spectrum is localized in time; and (3) 
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the system-characterizing quantities are determined from the impedance spectra using 
equivalent circuit fitting in a time-resolution-optimized form. 



BRIEF DESCRIPTION OF THE DRAWINGS 

5 

Figure 1 is a plot versus time of an example of a structured noise voltage 
signal U(t). 

Figure 2 is a power spectrum density plot of structured noise voltage and 
current signals U(f) and 1(f). 
10 Figure 3 is a plot illustrating examples of a sequence of impedance spectra as a 

function of the frequency f and the time x for a non-stationary system. 

Figure 4 is a schematic of an equivalent circuit. 

Figure 5 is a block diagram representation of one embodiment of a highly time 
resolved impedance spectroscopy system. 
15 Figure 6 is a table that shows examples of parameters for the structured noise 

signal and data processing for the described simulations in the "Simulations" column, 
and for the described measurements in the "fast model" and "Gigaohm model" 
columns. 

Figure 7 is a graph illustrating an example of the time course of the resistance 
20 R2 of the non-stationary system of Figure 4 (simulation). 

Figure 8 is a graph of an expanded portion of Figure 7 illustrating a transition 
in the value of R 2 . 



WO 01/27605 



PCTYUS00/26765 



8 

Figure 9 is a graph illustrating an example the time course of the values of the 
equivalent circuit elements of Figure 4 measured with highly time resolved 
impedance spectroscopy. 

Figure 10 (a) is a graph illustrating an example of the real part of the 
5 underlying impedance spectra for two points in the time courses in Figure 9. 

Figure 1 0 (b) is a graph illustrating an example of the imaginary part of the 
underlying impedance spectra for two points in the time courses in Figure 9. 

Figure 1 1 is a graph illustrating an example of the time course of the values of 
the equivalent circuit elements of Figure 4, measured with highly time resolved 
1 0 impedance spectroscopy. 

Figure 12 (a) is a graph illustrating an example of the real part of the 
underlying impedance spectrum for one point in the time courses in Figure 1 1 . 

Figure 12 (b) is a graph illustrating an example of the imaginary part of the 
underlying impedance spectrum for one point in the time courses in Figure 1 1 . 
1 5 Figure 1 3 is a block diagram illustrating one embodiment of the highly time 

resolved impedance spectroscopy. 

DETAILED DESCRIPTION 

20 A method and apparatus to measure impedance and impedance spectra in 

rapid sequence thereby providing repeated or continuous characterization of the 
electrical properties of the system under investigation. Highly time resolved 
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impedance spectroscopy is able to measure the dynamics of non-stationary systems 
with enhanced time resolution. 

Three aspects of highly time resolved impedance spectroscopy are: (1) an 
optimized, frequency rich a.c, or transient, voltage signal is used as the perturbation 
5 signal; (2) non-stationary time to frequency transformation algorithms are used when 
processing the measured time signals of the voltage U(t) and current I(t) to determine 
a sequence of impedance spectra, where each spectrum is localized in time; and (3) 
the system-characterizing quantities are determined from the impedance spectra using 
equivalent circuit fitting in a time-resolution-optimized form. 

10 

Optimized Frequency Rich Perturbation Signal 

The optimized frequency rich perturbation signal may be of many different 
formats. A frequency rich perturbation signal may contain any desired amount of 

15 contributing frequencies which are superimposed in one signal. For example, a 
frequency rich perturbation signal may be the superposition of a specific number of 
sinusoidal oscillations, or voltage jumps, pulses, and noise signals. The higher the 
signal to noise ratio of the voltage perturbation and the current response signal, the 
lower the time expansion of the database that is required to attain a certain measuring 

20 precision. 

Structured noise may be used as the perturbation signal for highly time 
resolved impedance spectroscopy. Structured noise is defined as the superposition of 
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a finite number of sinusoidal oscillations. Figure 1 is a plot versus time of an example 
of a structured noise voltage signal U(t) (18). Figure 2 is a power spectrum density 
plot of a structured noise voltage and current signals U(f) (20) and 1(f) (22). The 
structured noise illustrated in Figures 1 and 2, can be represented mathematically by 
5 equation (2): 

U(t)=YU 0l -sm(co { t+(p,) with 5</7<50 (2) 

/=i 

In equation (2), U 0 ,i is the amplitude of the i-th sinusoidal oscillation with 
angular frequency co\ and cpi is the phase. A favorable signal to noise ratio may be 

10 obtained when the perturbation signal U(t) is tuned, or optimized, to the system and 
the focus of measurement. Optimization of perturbation signals has been the subject 
of some research. See G.S. Popkirov and R.N. Schindler: "A New Impedance 
Spectrometer for the Investigation of Electrochemical Systems." Rev. Sci. Instrum. 63 
(11), 5366-5372 (1992).; G.S. Popkirov and R.N. Schindler: "Optimization of the 

15 Perturbation Signal for Electrochemical Impedance Spectroscopy in the Time 
Domain." Rev. Sci. Instrum. 64 (1 1), 31 1 1-31 15 (1993). All the above publications 
are incorporated by reference in their entirety herein. 

Numerous optimization steps are an aspect of highly time resolved impedance 
spectroscopy including: the frequency band in which the impedance is measured; the 

20 number and frequency, amplitude and phase of the individual contributing 
frequencies. 
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Optimization of the perturbation signal includes selection of the frequency 
band in which the impedance spectrum is measured. Therefore, the frequency band of 
the perturbation signal should be selected to include the maximum spectral range to 
be covered by the measurement. However, when evaluating the impedance spectra 

5 one can use just selected frequencies of the n frequencies in the perturbation signal. 
This aspect is discussed further below. 

Optimization of the perturbation signal may also include selecting a desired 
number (n) and frequency position (coi) of the n individual contributing frequencies. 
Any desired number and frequency positions of the individual frequencies can be 

10 used. This maximizes for the user the adaptability of the high time resolution 
impedance spectrometer to the system under investigation. 

In one embodiment, the power applied by a structured noise signal to the 
sample depends on the number of the frequencies in the noise signal. Because many 
systems exhibit linear behavior only with small amplitude or low power perturbation, 

15 it is desirable to keep the number of frequencies as low as possible. On the other 
hand, the reproduction of the characteristic system response requires a minimum 
number of frequencies in the impedance spectrum. Because in impedance 
spectroscopy the impedance is usually determined over several orders of magnitude of 
the frequency, a desired number of individual frequencies per frequency decade is 

20 selected, if special system requirements do not require deviating from this desired 
number. In one embodiment, five individual frequencies per frequency decade is a 
suitable number. 
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In one embodiment, the distribution of the frequency position of the individual 
frequencies is selected in a logarithmic uniform distribution over the entire frequency 
band. In another embodiment, the distribution of the frequency position of the 
individual frequencies selected is varied from a logarithmic uniform distribution in a 

5 manner to avoid the formation of harmonics. Preventing harmonics helps to prevent 
the corruption of the measured impedance spectra by nonlinear system responses or 
the excitation of higher harmonics. 

Another aspect of the perturbation signal that may be optimized is the 
amplitude U 0/i of the individual frequencies. Any desired amplitude of the 

10 individual frequencies may be selected. In one embodiment, to optimally exploit the 
linear range of a measuring amplifier, and attain an optimum signal to noise ratio for 
each frequency in the system, the amplitudes of the individual frequencies are adapted 
to the measuring situation. In one embodiment the amplitudes of the individual 
frequencies in the perturbation signal are such that the power of the individual 

1 5 frequencies in the perturbation signal is constant (good for use with strongly nonlinear 
systems). In another embodiment the power of the current response of the excited 
individual frequencies is constant (reduction of the influence of non-linearity of the 
measuring amplifier, see Figure. 2). In yet another embodiment, the power of the 
individual frequencies of the perturbation signal and the current response reveals 

20 minimal differences (optimum signal to noise ratio). 

Yet another aspect of the perturbation that may be optimized is the phase cp ± of 
the individual frequencies. Any desired phase angle of the individual frequencies 
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may be selected. In one embodiment, a linear system response maybe ensured by the 
lowest possible overall amplitude of the perturbation signal. The phases of the 
individual frequencies in the structured noise are selected so that the overall amplitude 
of the signal is minimal while the amplitudes of the individual frequencies remain the 
5 same. Hence the power of the individual frequencies is retained despite reduced 
overall amplitude, or the power of the individual frequencies is maximized and the 
signal to noise ratio is improved, as the overall amplitude remains fixed. 

Non-Stationary Time to Frequency Transformation Algorithms 

10 

An aspect of highly time resolved impedance spectroscopy is the 
determination of impedance spectra, which are localized in time, with reduced 
averaging time in high repetition rate. Non-stationary time to frequency 
transformation algorithms are used in highly time resolved impedance spectroscopy. 

1 5 The averaging time corresponds to the prolongation of the database that is required 
for determining an individual impedance spectrum. Sequences of spectra X(gd), which 
are localized in time, are determined from the time signals X(t) by applying window 
functions spanning a small period of time of the measuring signal X(t) (either U(t) or 
I(t)) and repeating this procedure, applying window functions shifted by brief, or 

20 short, intervals. The time signals, including the a.c. voltage signal U(t) applied to the 
sample and the current response I(t) of the sample, are continuously or partly 
continuously recorded over a given period only limited by the storage capacity of the 
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used data storage media. A partly continuously recorded data set may be defined as a 
repeated recording of finite data sets of any length with arbitrary, or selected, 
interruption intervals between consecutive data sets. 

Examples of non-stationary, time to frequency transformation procedures 
5 include sliding short time Fourier transform, wavelet transformation and Wigner or 
Wigner-Ville distribution. Two of these methods will be described further below: (1) 
sliding short time Fourier transform; and (2) wavelet transformation. 

In sliding short time Fourier transform the time to frequency transformation is 
carried out corresponding to the formula for continuous short time Fourier transform: 



In equation (3) x(t) corresponds to the measured time signal, g(t) is a window 
function whose characteristic is adjusted to optimize the spectrum X(t,co) that is 
obtained as a result of the transformation step. See L. Cohen: Time-frequency 
15 Analysis (Prentice Hall PTR, Englewood Cliffs, NJ, 1995) incorporated by reference 
in its entirety herein. Various possible window functions have been described. See 
F.J. Harris: "On the Use of Windows for Harmonic Analysis with the Discrete Fourier 
Transform." Proceedings of the IEEE 66 (1), 51-83 (1978) incorporated by reference 
in its entirety herein. 

20 The resulting impedance spectrum Z(t,co) is obtained according to equation (1) 

from the quotients of the spectra of the voltage signal U(x,co) and the current signal 



10 




(3) 
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I(x,co). This impedance spectrum is assigned to the time x, that is, it characterizes the 
investigated sample at time x of the measurement. The spectral information is hence 
localized at time x by considering, and weighting, only a section of the overall data 
record around time x for the time to frequency transformation. The impedance 
5 spectrum yields the complex impedance depending on the angular frequency co=2*7r*f 
where f is the frequency. A sequence of impedance spectra is obtained by repeatedly 
calculating impedance spectra according to the described procedure where the time x 
shifts by interval Ax for each of the subsequent spectra. 

Figure 3 is a plot illustrating examples of a sequence of the real (32) and 

10 imaginary (34) impedance spectra as a function of the frequency f and the time x for a 
non-stationary system. In the manner described above, the measured data records are 
completely analyzed when x is shifted as much as possible by interval Ax in the data 
record from the start of measurement. If the interval Ax is smaller than the 
prolongation of the window function g(t), one speaks of sliding short time Fourier 

1 5 transform since the window functions g(t-x) applied to the data records overlap within 
the sequence. The interval Ax in which the window function shifts between each 
transformation, defines the time resolution of the procedure. 

To attain a high time resolution, Ax must be correspondingly small. Likewise, 
the time over which g(t) is not zero must be as small as possible to reduce the time 

20 averaging in the time to frequency transformation. The window function g(t) defines 
the maximum peak-to-peak resolution At pp of short time Fourier transform provided 
by the expression: 
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(4) 



When the measured data records of U(t) and I(t) exist in discrete form, for 
5 example, in the case of computer-supported, or digitized, data acquisition, the discrete 
algorithm of the short time Fourier transform of equation (5) is used instead of 
continuous short time Fourier transform. 



10 number of observed data points, and N the total number of data points. All other 
quantities correspond to the definition in continuous short time Fourier transform. 
The shift interval Ax typically is a whole-number multiple of T a . The frequency 
resolution Af and the upper frequency limit f max in each impedance spectrum are 
subject to the sampling theorem expressed in equation (6). See E. Schriifer: 

15 S i gnal verarbeitung . (Hanser, Munich: 1990), incorporated by reference in its entirety 
herein. 



X d {T,co)=^{nT a )ginT a -x\e~^ ( 5 ) 



n=0 



T a is the sampling interval for data acquisition, n the running variable for the 



1 



and 



1 



max 



2-71 



(6) 



a 
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where N w is the number of data points above which the window function g(nT a ) is not 



In wavelet transformation the time to frequency transformation is carried out 
5 corresponding to the formula of continuous wavelet transformation: 



See O. Rioul and M. Vetterli: "Wavelet and Signal Processing" in Signal 
Processing Technology and Applications . Edited by J.G. Ackenhusen (The Institute 
10 of Electrical and Electronics Engineers, Inc., New York: 1995) p. 85-109 incorporated 
by reference in its entirety herein. 

Instead of the variable of the angular frequency co, the scale s is normally used 
for wavelet transformation. This is due to the fact that the pulse responses of the used 
wavelet functions h s (t) scales with s as e.g. is the case with a basic or prototype 
1 5 wavelet of the formula: 



zero. 




*,(0=-U<i) 




(8) 



s 



One possible form for the window function h(t) is a modulated window 
function g(t) as is used in short time Fourier transform: 
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h{t)= g(t)e~ 1 ^ , ( 9 ) 

In equation (9), co 0 corresponds to the modulation frequency of the basic, or 
prototype, wavelet. By scaling the basic function for the time to frequency 
5 transformation, the time resolution is not set for the entire spectrum but varies with 
the analyzed frequency or scale. This has an advantage that the time averaging is 
reduced by the time to frequency transformation at high frequencies by reducing the 
database. This increases the time localization and permits a substantially higher time 
resolution at this range. In addition, the considered database at low frequencies is 

10 expanded so that the contribution of these frequencies can be integrated into the 
spectral information. Overall, the information content of the determined sequence of 
impedance spectra may be greatly increased. The time resolution can be optimally 
adjusted by selecting the suitable shift interval At for a sequence of impedance 
spectra obtained by wavelet transformation, for either the entire frequency bandwidth 

1 5 or just selected frequencies of interest. 

In the case of discrete data records, a corresponding algorithm is used for 
discrete wavelet transformation or wavelet series expansion. 

The highly time resolved impedance spectroscopy offers many advantages 
over conventional impedance spectroscopy including the use of reduced impedance 

20 spectra. The impedance spectra obtained by short time Fourier transformation, e.g. 
intrinsically includes N w /2 + 1 frequencies which can be reduced to the n frequencies 
that are contributing to the structured noise. Because the discrete algorithms only 
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provide impedance values for discrete equidistant frequencies, frequencies are 
selected that are closest to the excited frequencies. This step eliminates most of the 
background noise of the measurement from the impedance spectra. Only the 
frequencies with the best signal to noise ratio are contributing to the impedance 
5 spectra, allowing the use of window functions which span a smaller period of time, 
resulting in a lower averaging time for the spectral information. 

Equivalent Circuit Fitting in a Time-Resolution-Optimized Form 

1 0 After determining the reduced impedance spectra, the spectra may be analyzed 

by adapting suitable equivalent circuit models to determine the system-characteristic 
quantities and their time course. The elements that compose the equivalent circuits are 
derived from physical models for the dominant processes occurring in the system and 
frequently correspond to the behavior of ideal electrical components such as resistors, 

1 5 capacitors and inductors. In the impedance spectrum, distinct processes dominate the 
impedance response of the system in different ranges of frequencies. In highly time 
resolved impedance spectroscopy, either the measured and reduced impedance spectra 
are evaluated completely or in certain ranges of the spectra. In principle, you can 
determine whether or not each frequency contributing to the reduced impedance 

20 spectrum is used for evaluation. 

An advantage of highly time resolved impedance spectroscopy is the 
possibility of multiple analyses of the same data records with different analysis 
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parameters, e.g. time resolutions Ax. For the multiple analyses of data records, first 
impedance spectra are determined with a longer window function leading to a long 
averaging time and a low time resolution, and therefore the reduced impedance 
spectra can be evaluated completely. The parameters of the stationary, i.e., constant 
5 processes are hence determined with a high degree of precision. 

The impedance spectra with a shorter window function, leading to a short 
averaging time and high time resolution, are determined in another analyses of the 
same data records. When the spectra are evaluated, only the areas of reduced spectra 
are used that characterize the dynamic non-stationary processes, and the parameters 
1 0 for the static or stationary processes in the equivalent circuit are set to the previously- 
determined values. 

The complex impedance values can be represented in different forms. For 
example, the impedance can be represented by complex coordinates or polar 
coordinates: 

15 

Z(co) = Re[z(©)] + / ■ Im[z(o>)] = Z 0 (©) -e l<p{(0 K ( 1 0 ) 

Likewise, the reciprocal impedance (termed admittance) or related quantities 
can be used such as the complex frequency-dependant dielectric constant in both 
coordinate forms. The significance of individual processes in certain frequency 
20 ranges of the impedance spectra varies with the coordinate form used. For analysis, 
highly time resolved impedance spectroscopy may use the electrical variables and 
coordinate form that significantly resolve the non-stationary processes in the highest 
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frequency range. The higher the frequencies at which the investigated processes can 
be analyzed allows for the selection of a smaller time period spanned by the window 
functions and a smaller shift interval Ax for use in the sliding short time Fourier 
transform. 

5 For example, a method for adapting the parameters of the equivalent circuit to 

the measured impedance spectra includes complex non linear least square fitting 
methods. A problem of the minimization algorithm for adapting the equivalent circuit 
to the impedance data (e.g. Levenberg-Marquardt, Powell's minimization method in 
multidimensions) is that the algorithm frequently finds local minimums instead of the 

10 global minimum in the fitting. See W.H. Press, S.A. Teukolsky, W.T. Vetterling et 
al Numerical Recipes in C (Cambridge University Press, New York, 1992) 
incorporated by reference in its entirety herein. 

To increase the probability of determining the global minimum in the fitting, 
different coordinate forms and variables can be combined. In highly time resolved 

15 impedance spectroscopy, sequences of impedance spectra may be evaluated by an 
automatic fitting routine or process. As a result, the time courses of the system- 
characterizing quantities are determined. 

The procedure described above is not limited to electrical impedance and can 
be used to determine mechanical impedance spectra such as in rheological measuring 

20 methods, magnetic and optical tweezers, quartz resonance balances and acoustic 
impedance measurements as well as other procedures in which the frequency rich 
perturbation of a system is related to the system response. 
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Exemplary Embodiment 

One embodiment will be described in detail to assist in understanding various 

5 aspects of highly time resolved impedance spectroscopy. Although one embodiment 
is described in detail, highly time resolved impedance spectroscopy may be embodied 
in other specific forms without departing from its spirit or essential characteristics. 
The embodiment describes a scenario that would facilitate the measurement, or 
determination, of the impedance of individual ion channels in supported lipid 

10 membranes, biomembranes. In another embodiment, the impedance of natural or 
artificial, freely-suspended or substrated supported lipid membranes with or without 
ion channels may be measured. The substrates bear planar thin-film microelectrodes 
that are in multiple electrode arrangements, facilitating parallel and sequential 
multiplex measurements on the substrate. The substrates can be fabricated from, for 

15 example, silicon substrates with metal or semiconductor electrodes. The substrates 
may be mounted in an electrochemical cell, of either one or multiple measuring 
chambers, in which liquid can be exchanged by a manual, or automated, liquid 
handling and control system with or without temperature control. Such a scenario 
corresponds to the measurement of an electrical two-terminal network that is 

20 simplified by the circuit depicted in Figure. 4. 

Figure 4 is a schematic of a typical equivalent circuit. The typical equivalent 
circuit of figure 4 has a resistor Ri (42) in series with the paralleled combination of 
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resistor R 2a (44), R 2 b (46) and C (48). A switch is in series with R 2a (44). The 
excitation voltage U(t) (50) is measured across the equivalent circuit, and the response 
current I(t) (52) flowing through the circuit is measured. Typical values of the 
components can be assumed to range from Ri ~ 100 kQ, R 2a , 2b ~ 1 GQ and C ~ 6 pF. 
5 Switching rates for the switch can be assumed to be e.g. 50 Hz. 

Figure 5 is a schematic representation of one embodiment. A structured noise 
of the a.c. voltage perturbation signal was calculated by a computer and sent to the 
memory of the multifunctional generator (54) (Analogic® 2030A). This device 
generates the perturbation signal and feeds it into the patch clamp amplifier (56) 

10 (HEKA® EPC8). 

When digital data records are used for the perturbation signal, the transitions 
between the discrete voltage steps should be small and smoothed because this 
reduces the level of background noise of the perturbation signal. In order to have 
small steps between two consecutive values of the perturbation signal, the sampling 

1 5 rate at which the data record is output by the function generator should be greater by a 
factor of about ten than the largest frequency of interest in the impedance spectra. 

To smooth the discrete steps, the perturbation signal can be filtered with a 
lowpass filter, with the corner frequency of the filter greater by a factor of about ten 
above the highest frequency of interest in the impedance spectra. The patch clamp 

20 amplifier (56) transmits the a.c. voltage perturbation signal via an external 
preamplifier (58) to the sample (60) under test. This preamplification unit 
simultaneously measures the actual a.c. voltage applied to the sample and the current 
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response of the sample. The current signal is filtered, amplified and undergoes 
current to voltage conversion. Both time signals U(t) and I(t) can be monitored as 
voltage signals at the corresponding outputs of the patch clamp amplifier (56). Both 
signals are acquired by two channels of the A/D converter board (National 

5 Instruments Lab-NB) in a measuring computer (Apple Macintosh Ilfx) and saved in 
data records on the computer data storage media. 

The software necessary to operate the setup includes two programs. The first 
program performs the data acquisition on the measuring computer with the A/D 
converter board. It includes the data analysis software Igor Pro® (WaveMetrics) with 

10 an expansion package for data acquisition (NIDAQ Tools). The second program 
performs the signal processing and data analysis of the measured data records of U(t) 
and I(t). It is a custom-written C++ program for Macintosh and Windows operated 
PCs. The implemented algorithm of the sliding short time Fourier transform is based 
on a modified Cooley-Tukey FFT algorithm. 

1 5 To fit the equivalent circuits to the impedance spectra, a modified minimizing 

algorithm was derived from the method of Powell in multidimensions. For the above- 
specified application, the logarithms of the real part and the negative imaginary part 
of the impedance spectra (complex coordinates) were used for the data fitting by the 
minimization procedure. See W.H. Press, S.A. Teukolsky, W.T. Vetterling et ah 

20 Numerical Recipes in C (Cambridge University Press, New York, 1992). 

The performance characteristics of the described embodiment were quantified 
in simulations and test measurements. For the simulations, the corresponding current 
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response signals of an ideal system were calculated by the computer for the structured 
noise perturbation signals, and the generated signals U(t) and I(t) were evaluated 
using the signal processing and data analysis procedure of highly time resolved 
impedance spectroscopy. Figure 6 is a table that shows examples of parameters for 

5 the structured noise signal (66) and data processing (68) in the "Simulations" column. 
The simulations demonstrate that a membrane/ion channel system with characteristic 
quantities Ri = 100 kQ, R 2 switching between distinct values in a range of 0.8-20 GQ 
and C = 6 pF can be measured correctly by using the logarithms of the complex 
coordinates in a minimization procedure even in a narrow frequency band of 0.3-20 

10 kHz with only six individual frequencies. In general, for the analysis of impedance 
spectra of membrane ion channel systems a frequency band of 1 00 Hz to 1 00kHz may 
be used. 

The equivalent circuit of the membrane/ion channel system is provided by 
Figure 4, where the combination of R 2a and R 2 b is considered as R 2 . 

1 5 In Figure 7 the time course (72) of the switchable parallel resistor R 2 is shown 

as it was assumed for the simulation (74) (continuous line) and as determined by 
highly time resolved impedance spectroscopy (76) (circle). Figure 8 shows one jump 
(82) in R 2 enlarged from Figure 7. The different curves shown result from using 
different window functions g(t) in the short time Fourier transform. Window 

20 functions of the following characteristics were used: Harming 1.0 (Cos.) (84), Kaiser- 
Bessel with a=2.0 (KB2) (86), and Kaiser-Bessel with a=3.5 (KB 3) (88). The 
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analysis of the time course reveals, that in the absence of any signal corruption, the 
time course of the system-characterizing quantities determined by the short time 
Fourier transform corresponds to the convolution of the real (ideal) time course of this 
quantity with the window function used. The time resolution attained was Ax = 0.67 
5 ms, and the rise time i r was between 2.0 and 3.1 ms depending on the window 
function. 

For the test measurements, the circuit in Figure 4 was assembled by using 
electronic components. Such a circuit forms a model sample, which mimics the 
behavior of individual ion channels in supported lipid membranes. These model 

10 samples were measured by highly time resolved impedance spectroscopy. 

In a first test measurement, a model was used with the values of R 2a , 2b = 680 
KQ and C = 1 nF (no additional element for Ri), and the switch was switched at a 
frequency of 50 Hz. Figure 6 shows the parameters of the structured noise signal and 
the data evaluation in the column, "fast model." In Figure 9 the time courses of Ri 

15 (92) , R 2 (94) and C (96) are shown as determined with highly time resolved 
impedance spectroscopy, as well as an ideal or absolute time course of R2 (98). 

Deviations from the absolute values of up to 20% result from the insufficiently 
compensated filter effects of the setup. The step-like changes of R 2 are reproduced 
well. The attained time resolution was 1.1. ms, and the rise time was 8 ms. In Figure 

20 10 (a) and (b) the real and negative imaginary part respectively of the underlying 
impedance spectra (112) and (104) for two points in the time courses of Figure 9 are 
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shown. Not each of the contributing frequencies of the structured noise signal was 
used for the evaluation but just the limited range of 0.4-4 kHz. 

In a second test measurement, a model was used with the values of R 2a ,2b = 1 
GQ and C = 3 pF (no additional element for Ri), and the switch was switched 
5 manually at irregular intervals. Figure 6 shows the parameters of the structured noise 
signal and the data evaluation in the column "Gigaohm Model." In Figure 1 1 the time 
courses of Ri (112), R 2 (114) and C (116) as determined by highly time 
resolved impedance spectroscopy are shown. Again, deviations from the absolute 
values of up to 20% result from the insufficiently compensated filter effects of the 

10 setup. The step-like changes of R 2 are clearly resolved also at these high impedances. 
The attained time resolution was 4.4 ms. In Figure 12 (a) and (b) the real and negative 
imaginary parts respectively of the underlying impedance spectrum of one point (122) 
in the time courses in Figure 1 1 is shown. Just five contributing frequencies in the 
impedance spectrum were sufficient to determine three independent, partially 

1 5 dynamic variables from the impedance spectra. This proves the significantly higher 
information content of measurements using highly time resolved impedance 
spectroscopy even in a very limited frequency band in comparison to measuring 
resistance at only one frequency. 

Figure 1 3 is a block diagram illustrating one embodiment of the highly time 

20 resolved impedance spectroscopy. As illustrated in Figure 13, the basic setup for 
highly time resolved impedance spectroscopy requires just a few components. A 
main processor (132) unit controls signal generation and data acquisition on two 
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channels. The data for the structured noise signal are generated by the processor unit 
according to a stored algorithm, stored in a memory device (134) and the voltage 
signal is applied to an output after a D to A converter (136). The measured voltage 
signal U(t) undergoes A to D conversion, at converter (138), at an input channel and 

5 is transferred to the memory. At the second input channel, the current signal I(t) is 
received by a signal conditioner (142) and undergoes current to voltage conversion at 
an A to D converter (144). The output of the converter (144) is transferred to 
memory. The measured data of the voltage and current signal are processed and 
analyzed with a routine stored in a memory device (146), and the results are saved on 

10 the data storage media (140). In another embodiment, the two memory devices (134) 
and (136) are the same memory device 

As described above, the use of the highly time resolved impedance 
spectroscopy has many advantages over conventional impedance spectroscopy. Some 
of the advantages include: substantial improvement in the time resolution of 

15 impedance spectroscopy; time resolution down to a fraction of each individual 
frequency in the impedance spectrum; continuous measurement with fast response 
times; real-time measurements; optimization of the measuring procedure by adapting 
the perturbation signal and the data processing and analysis to the investigated 
system; repeated analysis of the measured data records adapted to the focus of the 

20 evaluation; highly flexible technique; investigation of dynamic processes that 
previously could not be measured with impedance spectroscopy such as opening and 
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closing individual ion channels in natural or artificial, free or supported membranes; 

and highly favorable signal to noise ratio. 

The foregoing description details certain embodiments of the invention. It will 

be appreciated, however, that no matter how detailed the foregoing appears, the 
5 invention may be embodied in other specific forms without departing from its spirit or 

essential characteristics. The described embodiments are to be considered in all 

aspects only as illustrative and not restrictive and the scope of the invention is, 

therefore, indicated by the appended claims rather than by the foregoing description. 

All changes which come within the meaning and range of equivalency of the claims 
10 are to be embraced within their scope. 
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We Claim: 

1 . A method of highly time resolved impedance spectroscopy comprising: 
generating an ax. voltage perturbation signal and applying it to a system; 
recording the ax. voltage perturbation signal and a current signal produced in 

response to the perturbation signal; and 

processing the recorded signals using a non-stationary time to frequency 
transformation. 

2. A method as defined in Claim 1, wherein the non-stationary time to 
frequency transformation is a short time Fourier transform. 

3. A method as defined in Claim 1, wherein the non-stationary time to 
frequency transformation is a wavelet transformation. 

4. A method as defined in Claim 1 ? wherein the non-stationary time to 
frequency transformation is a Wigner distribution. 

5. A method as defined in Claim 1, wherein the recording of the a.c. 
perturbation signal is continuous. 
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6. A method as defined in Claim 1, wherein the recording of the a.c. 
perturbation signal is partly continuously recorded. 

7. A method as defined in Claim 1, wherein the a.c. perturbation signal is 
structured noise. 

8. A method as defined in Claim 1, further comprising a high- 
performance signal amplifiers such as patch clamp amplifiers. 

9. A method as defined in Claim 1, further comprising a analog-digital 
conversion of the measuring signals. 

10. A method as defined in Claim 1, further comprising data processing 
hardware and software for real-time analysis of the impedance spectra. 

11. A method as defined in Claim 1 ? further comprising detection of the 
conductivity of natural or artificial freely-suspended or substrate-supported lipid 
membranes with and without ion channels. 

12. A method as defined in Claim 1, further comprising the use of a 
frequency range of about 100 Hz to 100kHz and the logarithmic analysis of the real 
and imaginary part of the impedance. 
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13. A method as defined in Claim 1, further comprising detection of 
compounds in solutions with a conductivity changing effect in regard to 
biomembranes. 

14. A method as defined in Claim 1, further comprising using a measuring 
chambers with temperature control and liquid handling. 

15. A method as defined in Claim 1, further comprising use of a complex 
nonlinear least square method to fit parameters of an equivalent circuit to the 
measured impedance spectra. 

16. A method as defined in Claim 1, further comprising the use of 
automatic liquid handling systems. 

17. A method as defined in Claim 1, further comprising the use of 
multielectrode arrangements. 

18. A method as defined in Claim 17, further comprising the use of 
multiple measuring chamber systems with parallel data acquisition and analysis. 

19. A method as defined in Claim 17, further comprising the use of 
multiple measuring chamber systems with sequential data acquisition and analysis. 
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20. A method as defined in Claim 1 further comprising the integration in 
automated production and process flow. 

21. An apparatus for measuring a highly time resolved impedance spectra 
comprising: 

a memory device configured to store an algorithm for generation of a 
perturbation signal and instructions to process and evaluate measured data; 

at least one analog to digital converter configured to receive a signal 
representing a perturbation signal and to output a digital value representing the 
perturbation signal, and to receive a signal representing a system response to the 
perturbation signal and to output a digital value representing the system response; and 

a processor configured to receive the digital values output by the at 
least one analog to digital converters, and to receive the algorithm and instructions to 
process and evaluate measured data from the memory device, and to output a 
perturbation signal and an impedance spectra, wherein the instructions to process and 
evaluate measured data are non-stationary time to frequency transformation. 

22. The apparatus as defined in Claim 21, wherein the non-stationary time 
to frequency transformation is a short time Fourier transform. 

23. The apparatus as defined in Claim 21, wherein the non- stationary time 
to frequency transformation is a wavelet transform. 
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24. The apparatus as defined in Claim 21, wherein the non-stationary time 
to frequency transformation is a Wigner distribution. 

25. The apparatus as defined in Claim 21, wherein the perturbation signal 
is continuously recorded. 

26. The apparatus as defined in Claim 21, wherein the perturbation signal 
is partly continuously recorded. 

27. The apparatus as defined in Claim 21, wherein the perturbation signal 
is structured noise, 

28. The apparatus as defined in Claim 21, further comprising use of a 
complex nonlinear least square method to fit parameters of an equivalent circuit to the 
impedance spectra.. 

29. The apparatus as defined in Claim 21, further comprising the use of 
multielectrode arrangements. 

30. The apparatus as defined in Claim 21, further comprising use of 
automatic liquid handling systems. 

31. The apparatus as defined in Claim 21, further comprising use of a 
multiple measuring chamber systems with parallel data acquisition and analysis. 
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32. The apparatus as defined in Claim 21, further comprising use of a 
multiple measuring chamber systems with sequential data acquisition and analysis. 

33. The apparatus as defined in Claim 21, further comprising a measuring 
cell or chamber for detection of the conductivity of natural or artificial freely- 
suspended or substrate-suspended lipid membranes with or without ion channels. 

34. The apparatus as defined in Claim 21, further comprising a measuring 
cell or chamber for detection of compounds in solutions with a conductivity changing 
effect in regard to biomembranes. 

35. The apparatus as defined in Claim 21, further comprising the use of a 
measuring cell or chamber with temperature control and liquid handling. 
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